Capillary-like Fluctuations of a Solid-Liquid Interface in a Non-Cohesive Granular 

System 
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One of the most noticeable collective motion of non-cohesive granular matter is clustering un- 
der certain conditions. In particular, when a quasi-two-dimensional monolayer of mono-disperse 
non-cohesive particles is vertically vibrated, a solid-liquid-like transition occurs when the driving 
amplitude exceeds a critical value. Here, the physical mechanism underlying particle clustering relies 
on the strong interactions mediated by grain collisions, rather than on grain-grain cohesive forces. 
In average, the solid cluster resembles a drop, with a striking circular shape. We experimentally 
investigate the solid-liquid interface fluctuations, which are characterized through the static and 
dynamic correlation functions in the Fourier space. These fluctuations turn out to be well described 
by the capillary wave theory, which allows us to measure the grain-boundary surface tension and 
mobility once the granular "thermal" kinetic energy is determined. Furthermore, both quantities 
are consistent with a simple order of magnitude estimation considering the characteristic energy, 
length and time scales, which is very similar to what can be done for atomic systems. 

PACS numbers: 05.40.-a 45.70.-n, 68.08.-p 



Granular matter, composed by a collection of macro- 
scopic particles that interact via dissipative contacts, is 
fundamentally out of equilibrium. When energy is in- 
jected into an athermal granular system, constant grain 
rearrangements are induced and accompanied by friction 
and inelastic particle collisions. In particular, dense gran- 
ular materials exhibit a large variety of interesting phe- 
nomena, under static conditions or in a dynamical state. 
Notorious examples are granular piles, avalanches, seg- 
regation, pattern formation, granular phase coexistence 
and jamming [ll-Q . Since few decades, extensive research 
has been conducted on granular phase transitions. One 
case of interest is the solid-liquid coexistence during an 
avalanche, where the characterization of the two-phase 
interface allows the study of the granular rheology 
Another system of interest is vibrated granular layers 
that reveal Faraday waves jaatterns that are preceded by 
a solid-to-liquid transition [4, 5] . 

Dry granular systems are usually considered to have 
no surface tension. However, several recent studies show 
that non-cohesive or very weakly cohesive granular mate- 
rials develop phenomena driven by surface tension, which 
can be low but not zero. Some remarkable examples 
are the Rayleigh- Taylor-like instability in tapped pow- 
ders and the interfacial instabilities in falling granular 
streams, in air and vacuum Q. In other cases, granu- 
lar systems are shown to behave as a zero-surface-tension 
liquid, as for particle sheets (analog to "water bells") cre- 
ated by a granular jet impacting a target Q and finger- 
ing in a granular Hele-Shaw system [lO[. The point is 
to understand how capillary-like features can emerge out 
of collections of particles that are known to be almost 
or completely non-cohesive. In some cases, an hydrody- 
namic derivation conducting a zero-surface-tension limit 
succeed in describing the observations [l^l; in other 



situations, a low effective surface tension depending on 
the granular system dynamics is estimated, either con- 
sidering the flow of interstitial air Q or nano- Newton 
cohesion forces in the case of very low external forcing . 

Here, we experimentally investigate a granular system 
composed by non-cohesive mono-disperse grains that are 
confined in quasi- 2D geometry. This geometry allows 
the characterization of both individual grain motion and 
collective behavior. When such system is vertically vi- 
brated, it can present a transition from a completely flu- 
idized state to the coexistence of a liquid state with solid 
(ordered) clusters Ill4l3l|. It has been recently shown 
that such coexistence is triggered by a negative compress- 
ibility, as observed in a similar gas-liquid van der Waals 
granular transition Il4l.ll5i , and that density waves prop- 
agate in the system |12| . In our previous work 13| we 
show that depending on the vertical height and filling 
density the transition can be of either first or second 
order type. For both cases density fluctuations do not 
show strong variations at the transition, but local order 
varies strongly, either abruptly or continuously respec- 
tively, with a critical-like behavior in the second case. 

In this Letter, we stand in the stationary regime of 
phase coexistence in the case where the transition is con- 
tinuous. Above a critical driving amplitude a unique solid 
cluster is observed. We focus on the liquid-solid-like in- 
terface characterization. In average, the granular cluster 
exhibits a striking circular shape, like a drop, modulated 
by important grain-boundary fluctuations (see Fig. [T]). 
Our approach is similar to the one proposed in (3] , con- 
sisting in an analogy with condensed matter at molecu- 
lar scale where thermal agitation induces capillary waves 
that deform an interface. Thus, we aim to discuss to what 
extent it is valid to use some concepts of kinetic theory 
to describe our macroscopic athermal granular matter, 



2 




FIG. 1. (color online) (a) Typical image showing a solid-like 
cluster in coexistence with a liquid-like phase. The continuous 
line shows the grain-boundary detection, (b) Time averaged 
solid-liquid-like interface {R{6)) from about 7000 images. The 
continuous line shows the circle of radius Ro — 27.7 mm ob- 
tained from angle averaging {R{9)). 



in order to express the strong interactions mediated by 
grain collisions in terms of effective surface tension and 
grain-boundary mobility formulations. 

The experimental setup is the same one reported previ- 
ously It consists of a container confining a layer 
oi N = 11504 stainless steel spherical particles in a shal- 
low box with transverse dimensions = Ly = lOOd and 
height Lz = 1.94(i ± 0.02d, where d = 1 mm is the parti- 
cle diameter (configuration 2 of our previous study [l3l | ) . 
The top and bottom glass plates of the cell are directly 
in contact with the particles and the side walls. Particles 
are illuminated from below with an array of light emit- 
ting diodes and a high speed video camera enables par- 
ticle detection and tracking. The surface coverage is de- 
fined by the filling fraction </) = Nnd'^/AL'^ = 0.904. We 
submit the system to vertical sinusoidal vibrations, with 
displacement z{t) = ^sin(wt), where A is the vibration 
amplitude. Its motion is controlled with a piezoelectric 
accelerometer. Special care is devoted to the control of 
the horizontality. For this configuration, a second-order 
solid-liquid-like transition occurs when the dimcnsionlcss 
acceleration F = Au'^ /g reaches a critical value Fc ~ 5.1 
(/ = 2Tr /lu — 80 Hz). In the framework of the solid 
cluster's characterization, the present study is done at 
F = 6.30 ± 0.03 and / = 80 Hz. 

Figure [T] displays typical results of the interface de- 
tection. To distinguish the liquid phase from the solid 
phase, wc use a Voronoi cell's area criterion. For each 
particle, if its Voronoi area is lower than a critical value 
determined empirically [3], it is considered in the solid 
phase; otherwise, its in the liquid phase. Then, using 
coarse-grained polar coordinates {AO ~ 1°), we detect 
360 interfacial particles in each image. The origin of the 
polar coordinate system is fixed at the time-averaged cen- 
ter of mass of the solid-like particles. We observe that 
the time-averaged interface exhibits a circular shape with 
a mean radius of Ro/d — 27.7 ± 1.5. 



In this study, we consider a curve-driven interface be- 
havior for the steady state cluster. In condensed mat- 
ter, from the classical solid to liquid phases, any inter- 
face is microscopically rough due to the competition be- 
tween thermal energy and minimization of surface area 
[l7| . The grain-boundary evolution is involved in pro- 
cesses such as cry stallization in solution [3l or dendritic 
solidification [l9|, and controls structural and mechan- 
ical properties of many materials (l^l- The first inter- 
facial parameter studied is the grain-boundary stiffness, 
7 = 7 -|- 7", where 7 is the surface tension and 7" its 
second derivative with respect to the spatial coordinate. 
In our case, we make a small slope approximation for the 
interface, 7 7. The second interfacial parameter is the 
grain-boundary mobility M, defined by ^ = M-fK, with 
V the grain-boundary velocity and k the interface cur- 
vature. Inspired by theoretical and numerical studies on 
interfacial properties of molecular systems [21I [2^ . and 
on an experimental study of colloidal crystals [23|, we 
attempt to obtain these physical quantities applying a 
capillary wave description. 

In analogy to the capillary theory in condensed mat- 
ter, we will assume that there is a functional, analog to 
the free energy that is minimized in the stationary state 
and allows to obtain the dynamics close to the station- 
ary state. This assumption, although no fully justified in 
non-equilibrium systems, is made for simplicity and veri- 
fied a posteriori as its predictions are consistent with the 
experimental results. To follow the analogy with equi- 
librium systems, this functional will be refered as non- 
equilibrium free energy. 

First, we consider the interface contribution to the 
non-equilibrium free energy E^. In two dimensions, it 
is related to the cluster's arc length and by an effective 
surface tension 7 such as = ^ J^^ ^B? + {deRY d9, 
with R{9, t) the cluster's radius. Because of the system 
finite size, which implies a finite number of particles and 
a finite stationary radius, an additional mass term has to 
be added to obtain the total non-equilibrium free energy 
(for details see [l^l). The radius fluctuations are defined 
as 6R{6,t) = R{6,t) — i?o, where Rq is the mean radius 
in time and in space. Its Fourier representation is 

00 

dR{9,t)= 5R,n{t)exp{ime). (1) 

m—~oo 

For small radial fluctuations, it is direct to show that 
energy fluctuations obey 

SE = ^\6Ilo\' + ^ (|W?rp + (2) 

|m|>2 

where the first term corresponds to changes in size, the 
second to changes in position, which exists because the 
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FIG. 2. (color online) (a) The probability density function 
(PDF) of velocity fluctuations (open and solid symbols for 
radial and tangential components respectively): in the liquid- 
like domain (■), for the grain-boundary particles (•) and in 
the solid- like domain (a). The continuous line is a Gaussian 
fit to the PDF of the grain-boundary particle velocity fluctua- 
tions, (b) Average horizontal kinetic energy spectrum, {Km) 
versus m, in semi-log scale. The inset shows the linear plot 
for m = 0, .... 35. 



translational symmetry is not perfect in the experiment 
and the cluster has a tendency to remain in the center 
of the box, and the third term corresponds to changes 
of the cluster's shape. Only the last term is related to 
the surface tension 7 and the other two introduce new 
cocfEcients, v and A, that should be measured as well. 
From here, the static power spectrum is obtained 



(|<5i?o| 



^ (A-o)i?o 
' ttA ' 

(|^R|rn|^2p) = 



{K\m\^2)Ra 
TTj{m'^ — 1) 



(3) 
(4) 



where () denotes time average. The quantity {Km) = 
iiVpTOpdw^P -I- is the average horizontal kinetic 

energy per mode, where Np is the number of particles 
at the solid-liquid interface, mp = 4.45 ± 0.01 mg is the 
particle mass and and vlL are the interface particle's 
velocity Fourier components [2^. {K^) is the equivalent 
of the thermal energy kBT/2 per mode at equilibrium. 
Although related, {Km) is not equal to the usually de- 
fined granular temperature Tg ~ ^nipv'^, because the 
number of active modes is not 2Np. 

Figure [2^ displays the probability density function 
(PDF) of particle velocities in the liquid (1) and solid 
(s) phases, as well as for the particles at the bound- 
ary (b). The PDFs are almost Gaussian, with devia- 
tions for large fluctuations. They are centered on zero 
and we verify isotropy between the horizontal velocity 
components. The granular temperature Tg, obtained 
from the variance of the velocity distributions, is not 
uniform. Indeed, because of the collisional and dissipa- 
tivc nature of particle interactions, higher particle vol- 
ume fraction implies higher dissipation and lower Tg. 
As expected, we obtain > T^ > T|. In Fig. [Ja 
we present the average horizontal kinetic energy spec- 



trum, {Km) versus m. We observe energy equiparti- 
tion between m = 2 and m « 30, with an average 
value Kcq = 12.7 ± 0.2 nJ. For comparison, the bound- 
ary granular temperature Tg w 7.5 nJ. The two lowest 
modes have larger energy, {Kq) = 18.4 ± 11.3 nJ and 
{Ki) = 15.3 ± 7.2 nJ, and the energy components de- 
crease continuously to zero for m > 30. 

Figure [3^ shows a very satisfactory agreement for the 
experimental power spectrum ^ with a l/(m^ — 1) 
tendency, within a physically relevant range for the 
wavenumber m = 6 — 35, corresponding to wavelengths 
2TTRo/m « 5(i— 29d. Fitting Eqn. ([4]) to the experimen- 
tal data gives an effective surface tension 7 = 4.6±0.3 /iN 
in two-dimensions, and using Lz as the third dimen- 
sion, we get 73D = 'j/Lz = 2.4 ± 0.2 mN/m. Addition- 
ally, A and v can be measured through the components 
m ~ and to = 1 of {\SRm\'^) and {Km), obtaining 
A = 0.30 ± 0.28 mN, and u = 0.23 ± 0.18 mN. 

For comparison, our measured surface tension 73D is 
about 30 times smaller than pure water's surface tension, 
but it is much larger than the surface tension 0.1 ^N/m 
estimated for a freely falling dry granular material [8|. 
In fact, this latter work demonstrates that nano-Newton 
cohesive interaction forces, measured by AFM, are re- 
sponsible for the surface tension. In our case, the physi- 
cal mechanism is not originated by grain-grain cohesion, 
but from dissipative collisions between particles. Ac- 
tually, the effective surface tension of our system can 
be estimated as 7 ^ "^g/d ~ 7.5 /iN, which implies 
73D ^ Tg/{dLz) ~ 4 mN/m. This result is consistent 
with the scaling found in a recent numerical study on 
3D crystallization of hard spheres, where the fluid-solid 
surface tension is 730 ^ ksT/cP' [25j . 

From the definition of the grain-boundary mobility M, 
it is direct to show that each mode obeys a Langevin 
equation (see details in (il]): 



d5R„ 
dt 



-—SR„ 



Mf^m{t), 



(5) 



where t„i = RQ/[AI^{m'^ — 1)] for \m\ ^ 2, and tq = 
Rl/{M\) and r±i = Rl/{My) for to = and to = ±1 
respectively, are the relaxation times and rjmit) are the 
Fourier modes of the noise term ri{9,t) in real space, 
which is assumed to be delta correlated. From here, the 
following expression for the dynamic correlation function 
in Fourier space is derived: 



{6Rmit)SR*miO)) = {\6Rm\^)e-'/' 



(6) 



Figure [Sb displays a good agreement for the predicted 
exponential decay. The data for t > 0.2 s and for to > 18 
arc very noisy and not represented here. On one hand, 
the intercept of the exponential decay for t ~ should 
be equal to \og{{\6Rm\'^)), and so vary linearly with 
log(m^ — 1). Although it docs show this behavior, it is 
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log(m^ - 1) t (s) - 1 

FIG. 3. (color online) (a) Static analysis using « 7000 images. The average static power spectrum of the radius fluctuations 
{\5Rm\^), over the average horizontal kinetic energy spectrum {Km), as a function of m? — 1 in log-log scale, (b-c) Dynamic 

correlation function analysis, (b) Two-time correlation function {5Rm{t)5Rm(0)) for several m, which show exponential decay, 
(c) Decay time versus {rn? — 1). The continuous line shows the linear law fit. The error bars represent the dispersion over 
two successive experiments. 



systematically lower than log((|5i?m(0)|^}). When com- 
puting this last quantity, we observe that the 1/(to^ — 1) 
dependency is consistent. Fitting Eqn. (|4|) gives the same 
result for 7 than previously determined, 7 = 4.3±0.3 /iN. 
On the other hand, increases linearly with (m? — 1), 
as shown un Fig. [3j:. Putting numbers, we obtain the 
grain-boundary mobility AI = 5.0 ± 0.8 m'^J~^s~^. We 
remark that for an atomic system simulation [2l| . the 
mobihty in two-dimensions scales as M ~ l^/ {ksTtc): 
with fc^T the thermal agitation, and tc characteristic 
length and time respectively. In our case, considering 
ksT = Tg, Ic ^ d and tc ^ d/y/JSv^, yields the same 
order of magnitude for the mobility Af w 8 m'^J~^s~^. 

Another way to obtain the grain-boundary parameters 
A, v, and A/ is from a square-displacement analysis for 
TO = and TO = 1. Indeed, from Eq. ([5|), the mean 
square displacement for each mode is 

where Ai?o,i = (5i?o,i(i) — (5i?o,i(0). For short times 
t "C To,i, a diffusive behavior is expected for each mode, 
(|Ai?o,'i|') - 2£)o,ii, with 7^0,1 = M(Ko,i)/{7tRo), which 
is the analog of the fluctuation-dissipation relation used 
in atomic simulations for flat geometries [U . 

Figure |4] displays (|Ai?m p) — Cm as a function of time 
for TO = and to = 1. Here, Cm is a parameter obtained 
from the fit (|Ai?mp) = ^^[1 — exp(—i/rm)] -Kcm, which 
reflects the fact that the Langevin equation for each mode 
does not describe the first ballistic regime for very short 
times. The continuous lines represent the fit for each 
mode, showing a very good agreement. The predicted 
saturation is observed for long times, and the interme- 
diate diffusive regime is also manifested. From the fit- 
ted values of Am and r,„ we obtain A = 0.7 ± 0.4 mN, 
1/ = 0.4 ± 0.2 mN, A/o = 1.3 ± 0.8 m^J-^s^^, and 
Afi = 2.5 ± 1.2 m^J-is"^ Both mobilities should in 
principle be equal, and the same as measured for larger 
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FIG. 4. (color online) Mean square displacement analysis for 
m = and m — 1, in linear (a) and log- log scale (b). Each 
curve is the result of the average of two realizations. The solid 
lines correspond to fits (|Ai?„ip) = Am[^ — Gxp{—t/Tm)] + Cm■ 



m. The agreement between these values and those mea- 
sured from the static and dynamic correlation functions 
seems good enough considering the large errors in each 
quantity. 

To conclude, we have demonstrated that the liquid- 
solid-like interface in a quasi-2D vibrated granular sys- 
tem can be characterized by grain-boundary's parameters 
such as surface tension and interfacial mobility. Both 
quantities are consistent with a simple order of mag- 
nitude estimation considering the characteristic energy, 
length and time scales, which is very similar to what can 
be done for atomic systems. The scaling of the effec- 
tive surface tension with the granular temperature Tg 
suggests that the particles' kinetic energy plays the role 
of the cohesive energy that originates capillary-like phe- 
nomenon in molecular liquids. It would be interesting to 
relate the particle kinetic energy to a collisional pressure 
in each phase. By doing so, we could handle the surface 
tension concept by thinking in terms of pressure differ- 
ence, as defined by the hydrodynamic law of Laplace. 
This perspective would allow to couple a particle's scale 
study (e.g. by accounting for collisions, cross section 
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and contact duration determination) to a macroscopic 
description and would arise the question of particle pres- 
sure in granular media, well-known for homogeneous gas- 
fluidizcd bed (26| , but still open for high density granular 
flows. 
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